Spectral properties of statistical mechanics models 
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The full spectrum of transfer matrices of the general eight-vertex model on a square lattice is 
obtained by numerical diagonalization. The eigenvalue spacing distribution and the spectral rigidity 
are analyzed. In non-integrable regimes we have found eigenvalue repulsion as for the Gaussian 
orthogonal ensemble in random matrix theory. By contrast, in integrable regimes we have found 
eigenvalue independence leading to a Poissonian behavior, and, for some points, level clustering. 
These first examples from classical statistical mechanics suggest that the conjecture of integrability 
successfully applied to quantum spin systems also holds for classical systems. 



PACS numbers: 05.50.+q, 05.20.-y, 05.45+b 

Since the work of Wigner p]] random matrix theory 
(RMT) has been applied successfully in various domains 
of physics S. Recently several quantum spin Hamilto- 
nians have been investigated from this point of view. It 
has been found |||| that ID systems for which the Bethe 
ansatz applies have a level spacing distribution close to a 
Poissonian (exponential) distribution, P(s) = exp(— s), 
whereas if the Bethe ansatz does not apply, the level 
spacing distribution is described by the Wigner surmise 
for the Gaussian orthogonal ensemble (GOE), P(s) — 
exp(— tts 2 /A). Similar results have been found for 2D 
systems jS). This suggests that the GOE describes prop- 
erly some properties of the spectrum of complex quan- 
tum systems. In this letter we extend the RMT analysis 
from quantum spin systems to models of classical statis- 
tical mechanics. In particular, we look at possible conse- 
quences of integrability on the spectral properties of the 
transfer matrices. 

At first sight it seems natural to start with the Ising 
model in two dimensions without magnetic field as an 
example of an integrable model, and then to add a mag- 
netic field. It turns out that the spectrum of transfer 
matrices of the Ising model for finite size leads to numer- 
ical difficulties as explained below. We then have chosen 
the case of the general eight-vertex model on a square 
lattice (which contains the zero-field Ising model as a 
special case) 0. Moreover, it is known that zero-field 
eight- vertex transfer matrices commute with the Hamil- 
tonian of the anisotropic XYZ quantum spin chain for 
certain relations between the parameters of the two mod- 
els . We shall use the notation of Ref . j|] to designate 
the eight admissible vertices and their respective Boltz- 
mann weights a, a', 6, 6', c, c', d, and d! . We consider 
the row-to-row transfer matrices Tn to build iteratively 
a periodic rectangular lattice by adding rows of length N 
with periodic boundary conditions. Therefore the par- 
tition function of a periodic rectangle of n rows of N 
sites is Z n} M = Tr(Tjy-). Note that there can be different 
expressions for the matrix TV, but all these expressions 
have the same value of the trace of the nth power for any 
n, and therefore they can be deduced from each other by 



a similarity transformation; the spectrum is indeed an 
intrinsic quantity which does not depend on any partic- 
ular choice of the transfer matrix. To perform the usual 
statistical analysis of the spectrum we need to have real 
eigenvalues. However, in general this is not the case and 
we will restrict ourselves to cases where the transfer ma- 
trix is symmetric. It is well known that the eight-vertex 
model can be mapped onto an anisotropic Ising model 
on a square lattice with diagonal interactions and four 
spin interactions around each plaquette. We again use 
the notations of Ref. B and introduce the five coupling 
constants Jh, J v , J, J' and J". The transfer matrix of 
the spin model can be chosen symmetric if J = J'. In 
terms of the Boltzmann weights of the transfer matrix it 



requires only that c = d 



d' (c = d and d = d 1 



not a restriction). This condition is verified for models 
without electrical field (i.e. when a = a' and b = b') and 
also for models with a field. So we are able to build sym- 
metric transfer matrices for integrable cases with a = a', 
6 = 6', and c = d = d = d', as well as for non-integrable 
cases with o ^ a', d ^ 6', and c = d = d = d' . 

Before presenting our results, we briefly recall some 
features of the RMT analysis, which is a statistical anal- 
ysis of the eigenvalues of a given matrix regarded as an 
ordered set. Firstly, one has to sort the eigenvalues ac- 
cording to the symmetry of the corresponding eigenstate. 
In contrast to quantum spin systems, transfer matrices 
posses a priori no continuous symmetry (as e.g. the SU(2) 
spin symmetry), but only space symmetries. For row-to- 
row transfer matrices these are given by the automorphy 
group of a single row (and not of the full lattice). This 
group is the set of permutations g of sites such that g(i) 
and g(j) are neighbours if and only if i and j are neigh- 
bours. To each g acting on the set of vertices one can 
easily associate a linear operator g acting in the config- 
uration space. Obvioulsy g and XV commute. It is then 
possible to construct a set of projectors onto invariant 
subspaces of TV- This amounts to block-diagonalize TV. 
This is not only a useful way of lowering the size of the 
matrices to diagonalize, but also the manner to sort the 
eigenvalues. The automorphy group of the periodic ring 
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of length N is the dihedral group Dn generated by a 
translation and a reflection. Elemenentary group theo- 
retical analysis can be performed to build the (N/2 + 3 
if N is even, or (N — l)/2 + 2 if N is odd) projectors 
onto the invariant subspaces. The transfer matrix of the 
zero-field model is also invariant under the reversal of 
all arrows of the vertices and another projector has to 
be applied. Secondly, to find universal behavior within 
each invariant subspace, one needs to 'rescale' the eigen- 
values Ei in order to have a local density of eigenvalues 
equal to one. This operation is called the "unfolding" and 
produces the unfolded eigenvalues ej. The aim is to re- 
move the non-universal or system specific large scale vari- 
ations of the integrated density of states, and to study 
only the presumably universal short scale fluctuations. 
It amounts to compute an average integrated density of 
states N av (E) which is the smooth part of the actual inte- 
grated density of states. We then have e$ = N av (Ei). In 
the generic case, several methods can be used to compute 
N av (E): running average unfolding (local averaging of 
eigenvalues followed by interpolation), Gaussian unfold- 
ing (Gaussian broadening of each delta peak in the den- 
sity of states), and Fourrier unfolding (removal of short 
scale wave lengths using Fourrier transformation). 

The simplest quantity one calculates in RMT analysis 
is the distribution P(s) of the differences between two 
consecutive unfolded eigenvalues s< = e^+i — ej. For in- 
tegrable systems a Poissonian distribution is expected, 
since the ej are expected to be independent. In con- 
trast, for the simplest non-integrable systems the Wigner 
surmise is expected. Another quantity of interest is the 
spectral rigidity 0: 

/ , ra+L/2 \ 

A 3 (L) = (-min/ (N(e) - ae - bf de ) , 

\L a,b J a _ L/2 j a 

where (. . .) a denotes an average over a. This quan- 
tity measures the deviation from equal spacing. For a 
totally rigid spectrum, as that of the harmonic oscilla- 
tor, one has Ag sc (L) = 1/12, for an integrable (Pois- 
sonian) system one has Ag 01 ^) = L/15, while for the 
Gaussian Orthogonal Ensemble one has A3 E (L) = 
4y(log(L)- 0.0687) + <D(L- 1 ). It has been found that the 
spectral rigidity of quantum spin systems follows A3 01 (L) 
in the integrable case and A3° E (L) in the non-integrable 
case. However, in both cases, even though P(s) is in good 
agreement with RMT, deviations from RMT occur for 
A 3 (L) at some system dependent point L*. This stems 
from the fact that the rigidity A 3 (L) probes correlations 
beyond nearest neighbours in contrast to P{s). This is 
probably why the rigidity is much more sensitive to the 
parameters of the unfolding than the spacing distribu- 
tion. 

We have generated transfer matrices for different val- 
ues of the Boltzmann weights and linear size N = 16 
(resp. N — 14). This leads to matrix sizes of 65536 2 



(resp. 16384 2 ). Because of the eight-vertex condition, 
the transfer matrix couples only configurations with the 
numbers of up (or down) arrows having the same parity, 
so the matrix trivially seperates into two blocks. After 
projection the matrix splits up into 18 (resp. 16) sym- 
metry invariant blocks of which the largest has a size of 
2062 2 (resp. 594 2 ). We have tried several methods of 
unfolding. The results presented here are obtained us- 
ing either a Gaussian unfolding with a local broadening 
over five states, or a running average unfolding over ten 
states. In Fig. [fl we present a typical integrated density 




1e+06 2e+06 3e+06 4e+06 5e+06 6e+06 
E 



FIG. 1. The integrated density of eigenvalues N(E) (cir- 
cles) for the case a — 2.5, a' = 1.6, b = b' = 3, 
c = c — d — d! = l/v6 for a single symmetry invariant 
block. The insert shows how complex N(E) is even at a very 
fine scale. The full line is N a v{E). 

of eigenvalues together with the averaged curve. The 
spectrum is seen to be rather more complex than a usual 
spectrum of a quantum spin system, thus the unfold- 
ing has to be performed very carefully. Fig. || shows 
the probability distribution of the eigenvalue spacings 
averaged over all representations for two representative 
sets of Boltzmann weights. The diamonds correspond 
to a zero-field case with a = a' = 2, b = b 1 = 3, and 
c = c' = d = d! = l/v6- This point lies in an ordered 
region of the phase diagram, and the transfer matrix has 
eigenvectors of the Bethe ansatz form. The spacing dis- 
tribution is close to the Poissonian distribution. The plus 
signs on the same figure correspond to the case a = 2.5, 
a' = 1.6, b = b' = 3c = c' = d = d' = 1/V6. This 
point is in a region of parameter space where the Bethe 
ansatz does not apply. The spacing distribution is close 
to the Wigner surmise. In particular the level repulsion 
is clearly seen. To test further the idea that integrability 
leads to a Poisson law and non-integrability leads to the 
Wigner surmise we also have calculated the spectrum of 
transfer matrices for some particular points satisfying the 
free-fermion condition aa! + bb' = cc' + dd' || . In this 
case the results are less clear. In Fig. 13 we have chosen a 
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FIG. 2. The eigenvalue spacing distribution P(s) for two 
row-to-row transfer matrices of the eight-vertex model with 
N = 14. The diamonds correspond to a — a' — 2, b = b' = 3, 
c — c — d — d' = l/v6: this corresponds to an integrable 
point in the ordered region of the phase diagram without elec- 
trical field. The plus signs correspond to a — 2.5, a' = 1.6, 
b = b' = 3, c = c' = d = d' = l/v 7 ^ The full line is the 
Poissonian distribution while the dashed line is the Wigner 
surmise. 

free-fermion point within the generally non-integrable re- 
gion where an electrical field is present. The Boltzmann 
weights are a — 0.8, a' — 1/a, b = b' = \/2c 2 — 1, and 
c = c' = d = d! = 1/1.1. The spacing distribution is close 
neither to the Poissonian distribution nor to the Wigner 
surmise. In particular there is a peak at s = indicating 
level clustering. We have found the same phenomenon 
for the pure 2D Ising model which can be mapped onto 
a eight-vertex model satisfying the free-fermion condi- 
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FIG. 3. The eigenvalue spacing distribution P(s) for a 
point satisfying the free-fermion condition. The Boltzmann 
weights are a = 0.8, a' = 1/a, 6 = 6' = \/2c 2 — 1, and 
c = c' = d = d' = 1/1.1. Note the peak near s = and 
the Poissonian tail for s > 1. 

tion. This behavior is usually found together with a very 
involved density of eigenvalues and leads to numerical 



difficulties in the unfolding. A possible explanation of 
the peak at s = could be that the free-fermion model 
is a trivial non-generic model as for example the Hub- 
bard model at zero Coulomb repulsion. However, for 
some other values of the Boltzmann weights also obey- 
ing the free-fermion condition, the spacing distribution is 
much closer to the Poissonian distribution. This suggests 
as another possible explanation the existence of quasi- 
degeneracy leading to a Shnirclman peak at the origin 
H for this specific set of Boltzmann weights. This will 
be studied in detail in a forthcoming publication. From 
the above results we conjecture that the spacing distri- 
bution of eigenvalues of non-integrable models is close 
to the Wigner surmise corresponding to level repulsion, 
while for integrable models there is no level repulsion. In 
integrable systems there is level independence leading to 
a Poissonian spacing distribution, but with a tendency 
to level attraction in some cases. 

To go further in the RMT analysis, we present in Fig. || 
the spectral rigidity A3(L) for the same points in pa- 
rameter space corresponding to integrability and to non- 
integrability as in Fig. p[ The two limiting cases corre- 




FIG. 4. The spectral rigidity A3(L) for the same cases as 
in Fig. |. 

sponding to the Poissonian distributed eigenvalues (solid 
line) and to GOE distributed eigenvalues (dashed line) 
are also shown. For the integrable point the agreement 
between the numerical data and the expected rigidity is 
good up to a value L w 5. For larger values of L a 
saturation occurs showing the limitation of the model of 
independent eigenvalues. For the non-integrable case the 
departure of the rigidity from the expected behavior ap- 
pears at L ss 2, indicating that the RMT is only valid 
at short scales. Such behavior has already been seen in 
quantum spin systems We stress that these numeri- 
cal results depend much more on the unfolding than the 
results concerning the spacing distribution. 

In summary, we have numerically calculated the spec- 
trum of transfer matrices of the 2D eight-vertex model 
for various parameters. After having sorted and unfolded 



3 



the spectrum we have computed the eigenvalue spacing 
distribution and the spectral rigidity averaged over all 
representations. To our knowledge this is the first RMT 
analysis of transfer matrices of models in classical statis- 
tical mechanics. We have found that the non-integrable 
cases are well described at short scales by a the Gaussian 
orthogonal ensemble, while in the integrable cases the 
eigenvalues are mostly independent. We speculate that 
this is a general results. 
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